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Abstract 


The paper is devoted to the derivation of correction terms for the asymptotic formulae of the radiative 
transfer theory, which allow to increase their accuracy for smaller optical thicknesses (typically, down to the 
optical thickness equal to 5). Most of the examples are limited to the case of water clouds. However, the 
main results are easily generalized on other types of disperse media. 
© 2002 Published by Elsevier Science Ltd. 
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1. Introduction 


The geometrical thickness Z of many natural turbid media is much larger than the mean photon 
free path length / in the same medium. Clouds, snow, ice, foam, leaves and many other types of light 
scattering objects belong to this class. This is why the problem of light interaction with optically 
thick media attracted a lot of attention in the past. Convenient analytical equations for the reflection 
function of such media were obtained by Germogenova [1]. A specific case of weak absorption by 
optically thick media was considered by Rozenberg [2] and Zege et al. [3]. 

It was found [4] that the error of analytical solutions rapidly decreases with the optical thickness 
t= L/l, being smaller than 5% for t > 15 in the case of terrestrial clouds. The error rapidly grows 
with decrease of the optical thickness. This is due to the fact that in the asymptotic theory [1] terms 
of order t~” at n > 2 are neglected. The analytical derivation of next terms of the expansion of the 
reflection function on t is quite involved. 
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So the main task of this paper is to introduce empirical correction terms to analytical solutions, 
which broad the range of their applicability. Cases of nonabsorbing and absorbing media are con- 
sidered separately. 


2. Nonabsorbing media 


The reflection function of a nonabsorbing light scattering layer can be obtained from the following 
analytical equation [1,5-7]: 


RC, G 1, y) = R} (é, 1, Y) = tKo(č)Ko(n), (1) 


where € is the cosine of the light incidence angle vo, 7 is the cosine of the observation angle J, w 
is the relative azimuth, 


1 


= 2 
f a + 0.751(1 — g) c 
is the global transmittance, g is the asymmetry parameter and 
a=3 f KUE dG, 6) 
3 2n 1 
Kon) = = f ay f RAEE + MERE. (4) 
T Jo 0 


The function Ko(7) is called the escape function. It describes the angular distribution of photons, 
leaving the semi-infinite nonabsorbing scattering layer with light sources located at an infinite depth 
(the Milne problem). Note, that Eq. (1) can be easily modified to account for a Lambertian underlying 
surface with the albedo 4 [5]. We consider here, however, only the case of A = 0. 

We see that Eq. (1) reduces the calculation of the reflection function of a finite layer for large 
values of the optical thickness t to the calculation of the reflection function of a nonabsorb- 
ing semi-infinite layer R? (é, n, Y) with the same phase function as a finite layer under study. A 
user-friendly code for computing this function is freely available via Internet [8]. 

The approximate solution for R? (é, y, y) at nadir observation was obtained by Kokhanovsky [9]. 
It has the following simple form: 


0.37+1.94¢ | p(x — arccos č) 
ise AARO 
where p(@) is the phase function and @ is the scattering angle. It does not depend on the azimuth 


due to the symmetry of the problem. The approximate formula for the escape function is given by 
Sobolev [5]: 


Ko(€) = 5 (1 + 26), (6) 


which gives for «, using Eq. (3): « ~ 1.07. The accuracy of Eqs. (5) and (6) is better than 2% at 
¢ > 0.2. Thus, we have a chain of analytical equations (1), (2), (5) and (6), which can be used to 
calculate the reflection function of a finite optically thick layer at nadir observation. The accuracy of 


R? (eB (5) 
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this approximation is better than 5% at € > 0.2 and t > 15. To increase the accuracy of Eq. (1) for 
smaller values of the optical thickness, we found that an additional correction term in the expression 
for the transmission ¢ should be added. In particular, we derived, parametrizing the function (see 


Eq. (1)) 
R61) — RE 1) 


Č = ; 7 
KOKI) o 
where R(t, č, 1) is obtained with the exact radiative transfer equation solution, that 
=t to (8) 
where 
a — bë + cë? 
t = a — bẹ + cë“ (9) 


B 
and a= 4.86, b= 13.08, and c= 12.76. The value ¢* = t — tę should be used in Eq. (1) instead of t, 
if one likes to increase the accuracy of Eq. (1) for smaller t. The values of R? (é,1) and Ko(€) in 
Eq. (7) were obtained from Eqs. (5) and (6). Calculations were performed for a well-known Cloud 
C1 model of Dermendjian [10] at the wavelength 0.65 um, where light absorption by water droplets 
can be neglected. Note, that according to the reciprocity principle it follows from Eq. (9): 
2 bén + cn? éE 


c 


at nÆ£l. 
Finally, we have the following expression for the reflection function at nadir observation, which 


follows from Eqs. (1), (2), (5), (6), and (8): 


0.37 + 1.94 p(n —arccosč) 27 1 
Eee 4+) 49 \ 1.07 + 0.75t(1 — g) 


This equation reduces the calculation of the reflection function at nadir observation to that of the 
phase function calculation. Note, that the second term, containing the phase function, is smaller that 
the first one and can be often neglected (especially for rapid estimations). 

We study the accuracy of Eq. (11) in Figs. 1-3 for various values of optical thicknesses, effective 
radii of water droplets and incidence angles. The particle size distribution is given by 


f(a) = Ba® exp(—9a/aer), (12) 


where B is the normalization constant. The wavelength is equal to 0.65 um. 

In particular, the dependence of the reflection function on the solar angle is presented in Fig. 1 
for various values of the cloud optical thickness t and the effective radius of droplets aef. We see 
that Eq. (11) describes all important silent features of backscattered light. Namely, rainbow and 
glory phenomena present at solar angles around 40 and 0 degrees (scattering angles 140° and 180°, 
respectively). The rainbow is enhanced for the largest particles in Fig. 1. The minimum in glory 
scattering occurs for smaller scattering angles for smaller particles than it is for larger particles. The 
reflection function takes values from 0.1 till 1.1 in Fig. 1, while t increases from 5 to 100, which 


(10) 


T? 


R(t, é,1) = t) (1 +26) (11) 
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Fig. 1. The dependence of the reflection function of a plane-parallel homogeneous cloud layer with water droplets on the 
solar angle at nadir observation for different values of the cloud optical thickness, the effective radius of droplets and the 
wavelength 0.65 um (exact data—symbols, approximation (11)—lines). The droplet size distribution is given by Eq. (12). 
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Fig. 2. The error of approximation (11) as a function of the optical thickness for various solar angles and aef = 6 um. 
Other parameters as in Fig. 1. 


is a standard range of optical thicknesses for terrestrial water clouds. The approximation works well 
(see also Figs. 2 and 3) both for highly reflecting clouds and clouds with low reflection (R ~ 0.2). 

Note, that the dependence of the reflection function on the effective radius is more strong for 
smaller optical thicknesses, where single scattering still plays an important role. This dependence 
is almost disappeared for semi-infinite media due to highly developed multiple light scattering [9]. 
The glory scattering region is the most sensitive to the size of droplets. For instance, the reflection 
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Solar angle, degrees 


Fig. 3. The error of approximation (11) as a function of the solar angle for various values of the effective radius and 
t= 5,10. Other parameters as in Fig. 1. 


function can differ more by 50% for values of aef = 4 um as compared to the case of aef = 16 um 
(see the case of t=5) at small solar zenith angles. The difference, however, is smaller than 10% at 
t= 100 (see Fig. 1). 

Relative errors of Eq. (11) are given in Fig. 2 for different solar angles as the function of the 
optical thickness and dep = 6 um. They are smaller than 5%. The error is smallest at an incidence 
angle of 60°. The dependence of errors of the approximation on the size of droplets is presented 
in Fig. 3 as the function of the solar angle for optical thicknesses 5 and 10, where the largest 
error of the approximation is expected, and various values of the effective radius of water droplets, 
characteristic for water clouds, ranging from 4 to 16 um. Almost in all cases the error is smaller 
than 5%. It only weakly depends on the size of droplets. The errors are larger for Sun at the horizon 
(solar zenith angle is larger than 75°). However, this case cannot be treated in the framework of the 
plane—parallel geometry anyway. 

Summing up, the error of simple Eq. (11) is better than 5% at t > 5, Jo < 75° and aef =4-16 um 
for the case of water droplets. Thus, Eq. (11) can be used in a number of applications, including 
optical thickness determination from reflection function measurements. 

It should be pointed out that the reflection function of clouds is often used to their optical thickness 
determination [11], assuming a priori value of the effective radius of droplets. To study the influence 
of the deviation of the reflection function of clouds with droplet size distributions, having aef = 16 
and 4 um from that for the case of aef = 6 um, we have presented the differences 


_ R(aef = 6 um) ~ Ret = 16 ym) ee R(acf = 6 um) — R(der = 4 um) 


A 
R(aef = 6 pm) R(der = 6 pm) 


(13) 


in Fig. 4. Calculations were performed for a nadir observation. We see that these differences are 
in the range [—45%, 15%], depending on the size of droplets and the incidence angle. This can 
influence considerably the optical thickness retrieval. Thus, actual values of aef cannot be neglected 
in the cloud optical thickness retrieval algorithms. It follows from Eq. (11) that the dependence 
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Fig. 4. The differences 4ı(1,2), 42(3,4) (see Eqs. (13)) as functions of the solar angle. Other parameters as in Fig. 1. 
The values 4ı(2) and 4ı(4) correspond to the case of adjusted asymmetry parameters as described in the text. 


of the reflection function on the effective radius is determined by two separate factors, namely the 
phase function and the asymmetry parameter. To separate these dependencies, we present in Fig. 
4 also differences (13), assuming that the asymmetry parameter is correct and corresponds to the 
case of particle size distributions with aef = 16 and 4 um, respectively (see Eq. (11)), but the phase 
function corresponds to the case of droplets with aef = 6 um both for R(aef = 16 and 4 um). Then 
the differences are much smaller. They are within 5% for almost all solar angles (with exception of 
glory scattering, where it could reach 20%, see Fig. 4). This makes it possible to use the same phase 
function in the satellite cloud optical thickness retrieval. However, the value of g should be adjusted 
to the actual value of the effective radius of droplets, which can be obtained from measurements in 
the infra-red region of the electromagnetic spectrum. 


3. Weakly absorbing media 


The analytical result for the reflection function of weakly absorbing light scattering media was 
obtained by Rozenberg [2]. Namely, we have 


R(t, č., W) = Roœ(¢, n, Y) — T exp(—y — x)K (E)K (n), (14) 
where 
_ sh(y) 
~ sh(x + y) G2) 


is the global transmittance of a weakly absorbing layer and x = ,/3(1 — g)Bt, y= 4,/ 8/301 — g). 


The value of f gives the probability of absorption in a light scattering event. 
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It should be stressed that Eq. (14) should transform to Eq. (1) as 6 =0. This is not the case, 
however. We have from Eq. (15) at f =0: 
pen 1 
1 + 0.751 — g)t’ 
which differs from Eq. (2). 
Eq. (2) was derived from the asymptotical analysis of the radiative transfer equation. In derivation 
of Eq. (15) some phenomenological assumptions were made. This explains the difference. We propose 


to modify y in the denominator of Eq. (15) to make it consistent with Eq. (2). Then we have instead 
of Eq. (15) 


(16) 


_ shy) 
-shx + ay) 


Eqs. (17) transforms to Eq. (2) as 6 = 0. Moreover, we introduce the correction term (10) in 
Eq. (14). Then we have finally 


R(t, 6.11 W) = Rool Som W) ( 


(17) 


sh(y) exp(—x) 
sh(x + ay) 
where we also neglected the influence of the absorption on the escape functions. Eq. (18) reduces the 


calculation of the reflection function of a finite absorbing layer to that of a semi-infinite absorbing 
layer R.(é,4, Y). Note, that due to the relation [3] 


e) exp(— y)Ko(¢)Ko(1), (18) 


Rel, 1, Y) = R (č, N, W)exp(—uy), (19) 
where 
Ko(¢)Ko(1) 
= 20 
R(En ee 


the initial problem is reduced to the calculation of the reflection function of a semi-infinite nonab- 
sorbing layer. We found, using exact calculations, that the accuracy of Eq. (19) can be improved 
for larger f, modifying u in Eq. (19). Namely, one should use the parameter 


u* =(1—0.05y)u, (21) 


instead of u in Eq. (19). 

The accuracy of Eq. (18) in account with Eqs. (19)—(21), (5) and (6) is presented in Figs. 5—7 
for the case of the nadir observation at 2 = 1.55 um and different t and der. 

In particular, the cloud reflection function at the nadir observation is given in Fig. 5 for different 
cloud optical thicknesses and agg = 6 um. Both approximate results, obtained with Eq. (18), and 
the exact data are presented. The reflection function increases with the optical thickness, which 
corresponds to the dependence of the reflection function on the optical thickness is visible, where 
clouds do not absorb (see Fig. 1). 

Comparing Figs. 1 and 5, we see that the reflection function for absorbing clouds is smaller than 
for nonabsorbing clouds at t = 100. This is not the case for smaller optical thicknesses (see, e.g., 
cases of t=5 in Figs. 1 and 5). Such a behaviour can be understood taking into account that Figs. 1 
and 5 are prepared at different wavelengths. This results in different effective diffraction parameters, 
phase functions and asymmetry parameters for cases, presented in Figs. 1 and 5. In particular, the 
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Fig. 5. The dependence of the reflection function of a plane-parallel homogeneous cloud layer with water droplets on the 
solar angle at nadir observation for different values of the cloud optical thickness at the effective radius of droplets equal 
to 6 um, the wavelength equal to 1.55 um (exact data—symbols, approximation (18)—lines). The droplet size distribution 
is given by Eq. (12). 
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Fig. 6. The same as in Fig. 5 but for various effective radii and t = 5, 100. 


effective diffraction parameter of droplets 27der/2 is smaller at A= 1.55 um (see Fig. 5) as compared 
to the case of A2=0.65 um (see Fig. 1). Also the asymmetry parameter and, correspondingly, the light 
transmission at a given optical thickness and no absorption is also somewhat smaller at 1 = 1.55 um. 
This results in larger reflection for smaller particles, which is also confirmed by Fig. 1. Thus, almost 
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Fig. 7. The error of approximation (18), obtained using data in Fig. 5. 


the same values of the reflection function in Fig. 5 at t = 5 as compared to data in Fig. 1 at t=5 
are due to two factors, namely—the decrease of R due to absorption and increase of R due to the 
effect of the effective diffraction parameter of droplets. These two processes almost compensate each 
other in the case under study. 

The reflection function of a water cloud at nadir observation and various thicknesses of clouds 
and effective sizes of droplets is given in Fig. 6 as a function of a zenith solar angle. It follows 
from Fig. 6 that the reflection function of a cloud with the optical thickness of 100 and the size 
of droplets 16 um is close to the reflection function of the cloud with the optical thickness 5, but 
for the effective radius 4 um. This shows the importance the size of droplets for the reflection 
function calculation in the near-infrared as compared to the case of the visible part of the spectrum, 
considered in the previous section. It follows (see Fig. 7) that the error of Eq. (18) is smaller than 
5% at the wavelength 1.55 um and incidence angles smaller than 75°. Thus, this equation together 
with Eq. (11) can be used as a base of a semi-analytical algorithm for the size of particles and 
optical thickness retrieval from reflection function measurements at wavelengths 0.65 and 1.55 um. 

The complex refractive index of water at 1.55 um is 1.3109-0.00013487 and, therefore, the prob- 
ability of photon absorption f is quite small. It is equal to 0.0065 at aef = 6 um, for example (with 
g = 0.8214). It is interesting to study the influence of f on the error of the approximation (18). For 
this, we have calculated the phase functions (and also values of g and f, see Table 1) with the Mie 
theory, using the particle size distribution (12) at ag¢=6 um and the refractive index m= 1.3109 — ix, 
where y was taken equal to 0.00001, 0.00005, 0.0001, 0.0005, and 0.001. These calculated local 
optical characteristics then were used to find the reflection function both with the approximate theory 
and the exact radiative transfer code. We found (see Fig. 8) that for P < 0.0024 the error of the 
approximation is better than 5% for incidence angles smaller 75° and nadir observation. The error is 
equal to 5% also at B < 0.035(y < 1.2), but for smaller range of incidence angles 15° < Jo < 65°, 
which are quite typical for satellite observations. It amounts to approximately 15% for the case of 
the solar angle in nadir and 7=0.001. It should be pointed out that the value of y=0.0005 for water 
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Table 1 
The dependence of f,g, y on x 
x B g y 
0.00001 0.0005 0.8199 0.12 
0.00005 0.0034 0.8204 0.27 
0.00010 0.0047 0.8210 0.37 
0.00050 0.0226 0.8253 0.83 
0.00100 0.0435 0.8300 1.20 
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Fig. 8. The error of approximation (18) for various values of y = 0.12(1),0.27(2), 0.37(3), 0.83(4), and 1.2(5) (see also 
Table 1). 


corresponds approximately to the wavelength 2.312 um. This allows to use equations presented here 
for the case of water clouds both in visible and near-infrared. This range of wavelength is often 
used for the estimation of cloud properties from space. 


4. Conclusion 


We have proposed here modified equations for the calculation of the reflection function of optically 
thick media, which reduces the calculation of reflection function of weakly absorbing scattering 
layers to the calculation of the reflection function of a nonabsorbing semi-infinite layer. The code to 
calculate the reflection function of a semi-infinite layer is freely available via Internet [8]. At nadir 
observation the problem is reduced to the calculation of the phase function. Note, that the values 
of the phase function are usually small numbers in the backward hemisphere. Then the second term 
in Eq. (5) can be neglected as compared to the first one. This implies that we have an analytical 
equation for the reflection function at the nadir observation for incidence angles smaller than 75°. 
This analytical solution can be used, e.g., for the characterization of light scattering, absorption and 
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microphysical parameters of turbid media from spectral measurements of the reflection function. Our 
study was restricted to the case of water clouds. However, results with slight modification can be 
used for other light scattering media. The accuracy of analytical solutions obtained is better than 5% 
at solar angles less than 75° and nadir observation. 
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